Nucleon number dependence of longitudinal 
radii in ion-ion collisions as a signature of the 
onset of collective expansion 



Jan Pisut a ) 6 \ Neva Pisutova 6 ) and Petr Zavada c ) 



a)Laboratoire de Physique Corpusculaire,Universite Blaise Pascal, 
Clermont-Ferrand, F-63177 Aubiere,Cedex, France 
b) Department of Physics, Comenius University, SK-8^2 15 
Bratislava, Slovakia 
c) Institute of Physics, Czech Acad.Sci. Na Slovance 2, CZ-18040, Praha- 

Czech Republic 



Abstract 

In an attempt to disentangle the effects of nuclear geometry from 
those of expansion we study the dependence of the longitudinal ra- 
dius (Rl) on nuclear numbers of colliding ions within a simple model 
of multiple nucleon- nucleon collisions. The model has built in the nu- 
clear geometry but there is no dynamical expansion. In principle, what 
is significantly beyond models of this type can be considered as caused 
by collective expansion. Taking into account the present errors of data, 
oversimplified model and the absence of accurate information on Rl 
for pp interactions it is difficult to make strong statements about col- 
lisions of induced by light ions. On the other hand Pb-Pb data give 
Rl significantly higher than a geometrical picture could accomodate, 
presenting an evidence for the onset of collective expansion, or of onset 
of a new dynamical regime between O-Au and Pb-Pb collisions. The 
situation with S-Au and S-Pb collisions is less clear because of a rather 
large difference of data on Rl in these two cases. 



1 Introduction 



Heavy ion collisions provide the only way to laboratory studies of dense 
hadronic matter and hopefully also of theoretically predicted Quark- Gluon 
Plasma (QGP). Hanbury- Brown and Twiss (HBT) interferometry [1,2] 
brings information about dimensions of the homogeneity region [3,4,5] of 
particle production. 

Space-time evolution of proton-nucleus (pA) and of ion-ion (AB) collisions 
depends on details of dynamics and on values of important parameters, like 
the formation time of secondary hadrons,which are both only partially un- 
derstood. It is therefore impossible to predict accurately when the collective 
expansion of hadronic matter or of QGP sets on. In an ideal situation dy- 
namics of heavy-ion collisions would be known up to a few parameters which 
could be determined by HBT interferometry. Unfortunately the dynamics 
is unknown and because of that we shall start with simplest assumptions 
possible. 

We suppose that in the CERN SPS energy range collisions of protons 
and of light ions with nuclear targets consist of multiple nucleon-nucleon col- 
lisions,with dynamics close to that described in Ref.[6] . For earlier work in 
this direction, see Ref.[7] .In each of collisions,except for the last one, nucleon 
loses some fraction of its energy. In the last one it fragments roughly like in a 
pp collision. In those kinematical regions where the production of secondary 
particles is dominated by last collisions of individual nucleons,the model of 
Ref. [6] is close to predictions of the "wounded" nucleon [8] or "partici- 
pating" nucleon [9] models. The model of Ref. [6] admits a parton model in- 
terpretation, according to which fast partons, including valence quarks, do not 
participate in individual nucleon-nucleon collisions, what makes the picture 
consistent [10] with the experimental information on the A,B dependence of 
the Drell-Yan dilepton production. This model is not exactly that of mul- 
tiple nucleon-nucleon scattering in the whole rapidity domain. For secondary 
hadrons in the central rapidity region the picture might be close to that of 
independent nucleon- nucleon collisions, since slow hadrons in the nucleon- 
nucleon cms might be formed in between of successive collisions. 

The purpose of the present note is to discuss a formula for the dependence 
of Rl on nucleon numbers A and B of colliding ions, based on the picture of 
multiple nucleon-nucleon collisions. The formula might roughly correspond 
also to other models, provided that the collective expansion is rather small. 
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We shall then compare the formula with available data [11-14] of NA-44, 
NA-35 and WA-93 Collaborations obtained at the CERN-SPS and try to see 
whether the data can be understood only by nuclear geometry and fragmen- 
tation of nucleons or whether some dynamical expansion is required. We shall 
use two approaches. In the former one we assume that expanding hadronic 
matter is not present in proton and light ion collisions with nuclei in the 
CERN SPS energy region. The parameters in the model are then determined 
from this condition. Having parameters fixed in this way we compare the 
formula with data on Rl obtained in heavy ion collisions in particular in 
Pb + Pb and S + Pb ones. If values of Rl in these collisions are not larger 
than what follows from such an extrapolation of data from lower A's and B's 
then also S + Pb and Pb + Pb collisions would be most likely dominated just 
by nuclear geometry. If, on the other hand, longitudinal radii as measured in 
S+Pb and Pb+Pb are larger than such extrapolations, then a new dynamical 
mechanism,presumably caused by the collective expansion of matter, must 
be responsible for the increase oi Rl- 

In the latter approach we take the model literally and calculate the pa- 
rameters. Deviations from the model are supposed to give evidence for the 
presence of expansion. 

In both approaches we are trying to find a threshold at which Rl "de- 
couples" from simple geometrical dependence. Any point of decoupling or a 
point when one type of behaviour of Rl changes to another one is most 
likely related to the onset of a new dynamics of the collision. A threshold for 
" decoupling" from a simple geometrical dependence signals most probably a 
transition to a collective expansion of matter. 

Note that such thresholds should be correlated also with changes of pat- 
terns of other signatures, like enhancement of dilepton and photon production, 
increased J / ip suppression, appearance of azimutal asymmetries in noncentral 
collisions, etc. 
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2 Dependence of Rl on A and B in a model 
of multiple nucleon- nucleon collisions 

— * — * 

The correlation function for two identical pions with momenta ki and k 2 is 
written in the standard way 



C(fci,fc 2 ) = 1 + A 



j e ^-^)->(f; K)d 3 f 
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(1) 



where K = (ki + k 2 )/2 and p(f;K) is the density of sources of the cor- 
responding homogeneity region [3,4,5]. To simplify the discussion we have 
assumed that both pions have the same energy.In that case p(f; K) in Eq.(l) 
is integral over time of the space-time density distribution p{f , t; K) . The 

— * 

density of sources p(r; K) is asummed to be given, in some approximation, 

as 

p(r1 R) = exp (" 2ib " W~) (2) 

The Fourier transform corresponding to k\ — k 2 = q parallel to the beam 
(= z) axis becomes 

p(q) = exp{-^f) (3) 
and the correlation function in Eq.(l) is 

C(q) = l + \exp(-R 2 L q 2 ) (4) 
The longitudinal radius Rl is related to the mean squared value of z by 

R L 2 = (z 2 ) = J z 2 p(f; K)d 3 f (5) 

We shall further assume that in individual nucleon-nucleon collisions pions 
with small ^-component of momentum (in the nucleon-nucleon c.m.s.) are 
produced in a way similar to the one in collisions of free nucleons.The process 
goes pressumably via formation and decay of resonances in both cases. We 
thus assume that the density of sources for identical pions in a given nucleon- 
nucleon collision which took place in the point with coordinate z is given as 
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Here Ri{pp) is the longitudinal radius as obtained by HBT interferometry 
in pp collisions. The experimental information obtained by the AFS Collab- 
oration at the CERN ISR [15] indicates that R L (pp) is a bit larger than 
1/m.This result has been corroborated by further evidence obtained in less 
direct measurements [16] in studies of directional dependence of HBT inter- 
ferometry, based on the method proposed in Ref. [17]. The density of nucleon- 
nucleon collisions is parametrized as 

Pl {z) = j^eM-y (7) 

Where / is given by Eq.(6). The total density of sources is then given as 

p(z;K) = J p nn (z - z;K)pi(z - z)dz (8) 

where we have suppressed variables x, y.The Fourier transform of a convo- 
lution is the product of Fourier transforms and this leads via Eqs. (1,7,8,9) 
to 

C{q- K) = 1 + A exp (-(Z 2 (A B) + R 2 L {pp))q 2 ) (9) 

and finally to 

R 2 L (A,B) = R 2 L (p,p) + l 2 (A,B) (10) 
Where I 2 (A, B) is given in Eq.(7) and can be calculated as 

l 2 (A,B)=<z 2 > (11) 

The interpretation of Eq.(10) is simple. The experimentally observable longi- 
tudinal radius in AB collisions consists of two parts. The former, Ri(pp) is 
due to resonance decays or to equivalent dynamical reasons (production of 
clusters,etc.) and the latter is due to the geometry of the AB collision. 

3 Two approaches to the comparison of the 
formula with data 

We shall now compare the formula in Eq.(10) with data within the two 
approaches discussed in the Introduction. 
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The former approach is more qualitative. We assume that in A + B 
collision individual nucleon-nucleon interactions occur in a region with lon- 
gitudinal dimension 

I ~ crj- ((A - 1) 1/3 + (B - 1) 1/3 ) (12) 

Here r ~ 1.2/m,7 is the Lorentz contraction factor, c is a constant of the 
order 1 and replacements A — > (A — 1), .B — > (5 — 1) are unimportant for 
the data we shall discuss,but they facilitate the transition to pp collisions. 
The constant c may include some details of dynamics as well as the fact that 
the approximation ps 1.2tqA 1 / z underestimates [18] radii of light ions. 
Conventions used correspond to those of NA-44 and we have recalculated 
the NA-35 data accordingly. 

We shall consider the interference of identical pions with K z in the 
c.m.s. of nucleon-nucleon collisions. In this frame incoming nuclei appear 
contracted by the Lorentz factor 7 = s 1 / 2 /2M where s is the square of 
the nucleon-nucleon c.m.s. energy and M is the nucleon mass. For El = 
200AGeV,j ^ 10 and for E L = lQOAGeV we have 7^9. 

From Eqs.(10) and (12) we find 

R L \A, B) = R L \pp) + ((cr / 7 )((A - 1) 1/3 + {B - 1) 1/3 )) 2 (13) 

For A = B = 1 we obtain Rl(1, 1) = Rl{p-,p) as we should. 

The comparison of data on Rl(A, B) as obtained at the CERN SPS with 
Eq.(13) is presented in Fig. 1. The data corresponding to AB collisions were 
taken from the following sources: S + C from Ref.[ll], S + S from Ref.[ll]; 
p + Pb from Ref.[13]; S + Cu from Ref.fll]; S + Ag from Ref.fll]; O + Au 
from Ref.[ll]; S + Au from Ref.fll]; S + Pb from Ref.[13] and Pb + Pb from 
Refs.[13,14]. 

The straight line in Fig.l corresponds to the following values of parame- 
ters entering Eq.(ll) 

R L (pp) ~ 1.92/m; — w 0.36/m (14) 

7 

Since we have read the data of Ref.[13] from graphs and since the statis- 
tics is continuously increasing we do not give errors of the fit of data by 
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Eq.( 11), pointing out only that the determination of Rl{pp) by the data is 
less accurate than that of cr /7- 

The qualitative, and in our opinion the most important point, is easily 
visible in Fig.l. The data from S + C up to O + Au are compatible with the 
straight line whereas R^PbPb) is certainly much larger than what could be 
ascribed to the linear dependence obtained from proton- and light ion colli- 
sions with nuclear targets. We interpret this as an evidence for the presence 
of a threshold occuring between O + Au and Pb + Pb collisions at the CERN 
SPS energies. Fig.l contains also an indication of a possible discrepancy be- 
tween data on Rl on S + Au as given in Ref. [11] and on S + Pb presented 
in Ref. [13]. Difference between nucleon numbers of Au(197) and Pb(207) is 
hardly responsible for a large difference between the corresponding values 
of R L . If the true value of R L (SPb) is closer to Ref. [13] the threshold is 
probably around S + Pb collisions. 

The value of Rl{pp) ~ 1.92/m seems to be larger than expected, but since 
this is not very accurately determined by the "fit" we shall not try to analyze 
possible reasons for that. 

Taking 7 « 10 and r w 1.2 fm we find c » 3 what is also somewhat 
large. The value of c might be also influenced by the underestimate of radii 
of light ions by the expression Ra = roA 1 / 3 ,see Ref. [18]. 

Note that in the region where data in Fig.l are presented the dependence 
as given by Eq.(ll) is not much different from a linear relationship between 
R L (AB) and (A - I) 1 / 3 + (B — I) 1 / 3 . 

Arguments given above are admittedly rather crude. In spite of simplifica- 
tions used we shall make an attempt to interpret the results. The experimental 
value [13,14] of the square of longitudinal radius for Pb+Pb collisions is about 
40 fm 2 , whereas the extrapolation from proton- and light ion nuclear colli- 
sions gives about 20/m 2 .The expansion should be responsible for the differ- 
ence. Putting R 2 L ,exper ~ R 2 L,geom + R 2 expan We get R eX pan ~ 6/m.ShlCe (z 2 ) = 

Rl 2 we have (z 2 ) 1 ^ 2 expan ~ 5. 6fm. Assuming that for pions with K z pz the 
homogeneity region corresponds roughly to rapidity region —1/2 < y < 1/2 
than the time of expansion is roughly t expan pz ( z2 ) l JxpaJ v iv = V^) where 
v(y = 1/2) = tanh(l/2) is the velocity corresponding to y — 1/2. In this way 
we get the estimate t expan pa 11.2 fm/c. 

The latter approach is more rigid. Values of l 2 (A+B) entering Eqs.(10) 
and (11) are obtained from nuclear geometry without any free constants. We 
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have calculated values of I 2 (A + B) for different colliding nuclei by using 
the model of Ref.[19].In this model distribution of nucleons in nuclei is given 
by a standard Woods- Saxon density. The resulting values of I 2 (A + B) for 
central collisions in the c.m.s of nucleon- nucleon interactions and taking into 
account Lorentz contraction are as follows (in fm 2 ) 

l 2 (S + C) = 0.0U 1 2 {S + S) = 0.018 1 2 {S + Cu) = 0.022 

l 2 (S + Ag) = 0.029 l 2 [0 + Au) = 0.037 l 2 (S + Au) = 0.041 (15) 

l 2 (p + Pb) = 0.038 l 2 (S + Pb) = 0.042 Z 2 (P6 + Pb) = 0.055 

In these calculations we have not taken into account deceleration of nucleons 
by collisions. 

Very similar values of I 2 (A + B) are obtained in a simple model in which 
nucleus with radius R is replaced by a cylinder of length 2L with L = 
^/3/5-R.This relation is given by the requirement that < z 2 > is the same in 
both cases. 

Values in Eq.(15) are rather small with respect to R 2 i{pp) and the corre- 
sponding curve in Fig.l is practically flat. To draw the curve we need the value 
of R 2 L {pp). The interpretation of data is thus made difficult by the absence 
of experimental information on i?|(pp).But even admitting a very large, and 
probably unrealistic, value of R 2 L (pp), say up to about 4fm 2 the data would 
indicate a presence of expansion even in collisions like p+Pb. Assuming a par- 
ticular value of R\{pp) we can roughly estimate the expansion times as above. 
For Rl(pp) = Ifm and Ri{Pb + Pb) ~ 6fm we obtain t expan ~ 12.6/m/c for 
Pb-Pb interactions and with Rl(S + Pb) 4/m we have t expan 8.4/m/c 
for S-Pb collisions. 

4 Comments and conclusions 

Data on Rl obtained in ion-ion collisions at the CERN-SPS [11,12,13] give 
a valuable information on the space- time evolution of these collisions. We 
have tried above to analyze these data by using a very simple picture of ion- 
ion collision as a sum of nucleon- nucleon interactions. The picture leads to 
a simple formula Eq.(10).We have used this formula in two approaches. In 
the former we have introduced Rl{pp) and the constant c in Eq.(12) as free 
parameters. In this way data with incident light ions can be interpreted as 
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due to nuclear geometry, but Rl for Pb-Pb collisions is definitely higher than 
a geometrical model can accomodate what gives evidence for the presence 
of collective expansion of matter in Pb-Pb collisions. The case of S-Pb and 
S-Au collisions is difficult to interpret since there is an unexpectedly large 
difference between the two data. It turns out also that an accurate information 
on Ri{pp) is of primary importance. 

In the latter approach we have calculated the contribution of nuclear 
geometry to Rl(pp) i n a more rigid wayleaving Rl(pp) as the only unknown 
parameter. The resulting Rl(AB) as given by geometry is almost independent 
of A, B and practically equal to Rl(pp) .The value of Rl{pp) is thus a key to 
interpreting data on Rl for light- ion induced interactions. Pb-Pb data are 
again much higher than geometry can explain. Our estimates of the expansion 
time in Pb-Pb collision give t expan ~ 9/m/c. 

Study of the dependence of Rl on A,B and on incident energy can, in 
our opinion, reveal the changing dynamics of ion- ion collisions. This con- 
cerns in particular presence of "thresholds" corresponding to changes of the 
dynamical regime. 

The absence of accurate information on Rl(pp) and rather large errors 
of data make it difficult to see whether data indicate presence of one or two 
"thresholds". The first one might be rather low - somewhere between pp and 
S-C collisions,but this question depends crucially on the value of Rl{pp), and 
the second one between O-Au and Pb-Pb. The existence of the latter seems 
very probable. 
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Figure Caption The dependence of y = Rl 2 (A, B) as given in Refs.[ll- 
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14] on x = ((A— l) 1 / 3 + (B — l) 1 / 3 ) 2 . Straight line fit corresponds toy = a + bx 
with a = 3fm 2 and b = 0.13/m 2 .In making the fit we have excluded data 
on S + Pb and Pb + Pb collisions. Note the Pb + Pb point in the right upper 
corner. For discussion of the second approach to data the reader is invited to 
draw a line parallel with the x-axis and corresponding to his preferred value 
of R L {pp)- 
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